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A  Novel  Empirical  Orthogonal  Function  (EOF)-Based 
Methodology  to  Study  the  Internal  Wave 
Effects  on  Acoustic  Propagation 

GaSlle  Casagrande,  Yann  Stephan,  Alex  C.  Warn  Vamas,  and  Thomas  Folegot 


Abstract — ^Thi.s  paper  presents  a  novel  approach  to  synthesize 
realistic  environment  for  ocean-acoustic  parametric  studies.  In  its 
current  form,  this  methodology  applies  to  internal  waves  and  tides. 
Empirical  orthogonal  function  (EOF)  decomposition  is  applied  to 
a  temporal  series  of  temperature  profiles.  It  can  be  observed  that 
the  first  two  time-dependent  expansion  coefficients  are  dynami¬ 
cally  linked.  When  they  arc  plotted  one  versus  another  in  a  scatter 
diagram,  the  cloud  of  points  consists  of  a  crescent  shape  that  can 
easily  be  represented  by  a  polynomial  fit.  If  the  first  two  expansion 
coefficients  capture  enough  variability  in  the  temperature  profiles, 
the  EOF  modes  plus  the  polynomial  can  be  used  to  reconstruct  tem¬ 
perature  profiles  independently  from  the  set  of  data.  This  realistic 
synthesized  environment  can  then  be  input  to  acoustic  propaga¬ 
tion  models.  This  approach  is  applied  to  the  case  of  the  Messina 
Strait  in  which  internal  waves  are  known  to  be  intensive.  From  a 
short-term  series  of  temperature  profiles  collected  on  a  thermistor 
string,  range-dependent  profiles  along  and  across  the  strait  are  re¬ 
constructed.  The  acoustical  impact  study  is  conducted  with  the 
range-dependent  acoustic  model  (RAM)  parabolic  equation  (PE) 
model.  The  methodology  presented  in  this  paper  is  simple  to  run 
and  requires  a  very  affordable  set  of  data.  It  could  be  used  as  an  ef¬ 
ficient  alternative  to  ocean  and  acoustic  model  coupling  for  process 
studies  or  for  regional  studies  especially  in  poorly  known  areas  or 
highly  variable  areas,  where  it  is  difficult  to  obtain  good  sound- 
speed  profile  prediction  from  ocean  models. 

Index  Terms — Empirical  orthogonal  function  (EOF),  internal 
waves,  Messina  Strait,  ocean-acoustics  modeling. 


1.  Introduction 

INTERNAL  WAVES  (IW)  are  known  to  have  a  significant 
impact  on  the  acoustic  propagation  (e.g.,  [I]  and  [2]).  From 
a  geometrical  point  of  view,  the  vertical  displacements  of  the 
thermocline  can  modify  drastically  the  surface  duct  and  can 
also  locally  modify  the  conditions  of  propagation  if  the  acoustic 
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source  or  receivers  are  placed  in  this  Critical  layer  [3].  From  an 
energetic  jx)int  of  view,  the  acoustic  propagation  through  IW 
can  also  induce  focusing  effects  [4].  Anticipating  the  effects  of 
internal  waves  in  a  region  of  interest,  whether  scientific  or  mili¬ 
tary,  is  then  highly  interesting  for  tomdgraphers  and  oceanogra¬ 
phers,  who  may  want  to  design  and  optimize  inversion  schemes 
or  for  sonar  operators  and  decision  rpakers  who  may  want  to 
optimize  antisubmarine  warfare  (AS\V)  tactics  and  behaviors 
in  IW  fields. 

The  impact  of  IW  can  be  numerically  studied  coupling  IW 
models  and  ocean-acoustic  propagation  models.  Such  studies 
can  be  separated  into  two  categories^  process  studies  and  re¬ 
gional  studies.  As  far  as  process  studies  are  concerned,  the  pur¬ 
pose  is  generally  to  conduct  a  parametric  study  of  the  impact  of 
IW  in  a  set  of  acoustic  scenarios.  For  that  purpose,  simplified 
IW  models  can  be  used  as  feature  models,  relying  on  a  math¬ 
ematical  description  of  the  IW  or  process  models,  numerically 
resolving  the  underlying  physics  of  the  generation  and  propa¬ 
gation  of  IW,  whereas  in  regional  studies,  the  aim  is  to  provide 
a  realistic  environment  in  which  IW  i$  only  a  contribution.  For 
that  purpose,  regional  general  circulation  models  can  be  used 
provided  they  can  account  for  IW  predictions. 

We  propose  in  this  paper  a  novel  approach  which  can  be  con¬ 
sidered  as  a  compromise  approach  between  process  modeling 
and  regional  modeling.  The  motivation  of  this  work  is  indeed 
to  have  a  simple  efficient  way  to  build  realistic  environments 
when  one  wants  to  focus  on  a  particular  phenomenon  in  a  par¬ 
ticular  region.  Our  approach  is  data  driven  but  does  not  re¬ 
quire  a  full  knowledge  of  the  region  of  interest.  On  the  other 
hand,  our  approach  is  also  parametric  but  does  not  require  a  too 
complicated  process  model  or  a  too  pimple  parametric  model. 
The  method  proposed  in  this  paper  relies  on  the  decomposition 
of  a  simple  set  of  data  (in  practical,  a  short-term  temporal  se¬ 
ries  of  temperature  profiles)  into  empirical  orthogonal  functions 
(EOFs).  It  has  been  shown  in  [5]  that  in  the  case  of  I W  in  the 
Messina  Strait,  the  time-varying  expansion  coefficients  of  the 
first  two  EOFs  represent  over  99%  of  the  variance  and  are  dy¬ 
namically  linked.  When  plotting  the  expansion  coefficients  one 
versus  the  other,  the  scatter  diagram  has  a  well-defined  cres¬ 
cent  shape.  This  feature  is  confirmed  by  EOF  analysis  of  numer¬ 
ical  outputs  from  the  Lamb  model  [6].  The  signature  of  internal 
waves  in  EOF  expansion  coefficient  Scatter  diagrams  can  then 
be  used  to  detect  the  presence  of  internal  waves  in  regions  of 
the  ocean.  A  curve  can  be  fitted  to  the  crescent- shaped  scatter 
diagram  in  presence  of  internal  waves.  Such  a  curve  could  be  de¬ 
rived  for  various  regions.  This  led  to  building  a  representation 
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Fig.  I .  Bathymetry  of  the  Strait  of  Messina;  the  color  bar  indicates  the  depths  in 
meters.  The  blue  dot  represents  the  position  of  the  thermistor  string.  The  dotted 
lines  correspond  to  the  axis  of  simulations  in  the  two  acoustic  configurations. 
For  the  transverse  configuration,  the  source  is  located  close  to  the  Sicilian  coast 
and  emits  towards  the  Italian  coast.  For  the  axial  configuration,  the  source  was 
positioned  cither  in  front  of  Messina  and  emitting  southwards  (downslope  con¬ 
figuration),  or  at  the  south  part  of  the  strait  and  emitting  towards  the  sill  (upslope 
configuration). 


of  internal  wave  events  with  no  more  that  a  statistically  perti¬ 
nent  set  of  data.  A  parametric  study  enables  to  reconstruct  a  set 
of  plausible  range-dependent  environments  which  account  for 
the  actual  variability  in  the  considered  region.  By  applying  an 
acoustic  propagation  loss  model,  it  is  then  possible  to  antici¬ 
pate  the  acoustical  impact.  This  novel  methodology  is  applied 
to  a  data  set  from  the  COACH ’06  experiment  jointly  conducted 
in  November  2006  by  the  NATO  Undersea  Research  Centre 
(NURC)  and  the  Service  Hydrographique  et  Oc6anographique 
de  la  Marine  (SHOM). 

In  the  first  part  of  this  paper,  we  recall  the  data  set  and  de¬ 
scribe  the  variability  in  terms  of  EOF  decomposition.  We  then 
describe  the  method  to  reconstruct  the  synthetic  environments. 
To  illustrate  our  approach,  the  third  part  presents  a  parametric 
acoustical  study  of  the  effect  of  IW  in  the  Strait  of  Messina 
using  the  variability  obtained  by  EOF  decomposition.  Conclu¬ 
sions  and  perspectives  are  drawn  in  the  last  part  of  the  paper. 

11.  Data  AND  Statistical  ANALYSIS 

A.  Preliminary  Consideration  on  the  Strait  of  Messina 

The  Strait  of  Messina  is  a  rich  dynamic  area  where  many 
ocean  features  can  be  encountered.  It  separates  the  Italian  Penin¬ 
sula  from  the  island  of  Sicily  and  is  a  natural  connection  be¬ 
tween  the  Tyrrhenian  Sea  in  the  north  and  the  Ionian  Sea  in  the 
south.  The  strait  is  a  narrow  channel,  whose  smallest  cross-sec¬ 
tional  area  is  0.3  km^  in  the  sill  region.  There,  the  mean  water 
depth  is  80  m  (Fig.  1 ).  The  water  depth  increases  gently  in  the 
northern  part  and  rapidly  in  the  southern  basin.  Throughout  the 
year,  two  water  masses  are  encountered  in  the  strait:  the  Tyrrhe¬ 
nian  Surface  Water  (TSW)  and  the  colder  and  saltier  Levantine 
Intermediate  Water  (LI  W). 


Although  tidal  displacements  are  Very  small  in  the  strait  of 
Messina  (the  order  of  10  cm),  the  tidal  phase  opposition  be¬ 
tween  the  two  basins  can  lead  to  a  sea  surface  slope  of  I  -2  cm 
km*"^  thus  generating  very  strong  cuifents  as  high  as  3  m  s”^  in 
the  sill  region  during  spring  tides  [7],  [8].  The  strong  barotropic 
tidal  flow  over  a  steep  bathymetry  in  $  stably  stratified  environ¬ 
ment  represents  the  ideal  recipe  for  internal  wave  generation  [9]. 
Baines  [10],  [1 1]  has  demonstrated  the  generation  process  ana¬ 
lytically.  The  hydrodynamic  modeling  of  tidal  flow  over  steep 
topography  shows  that  the  tidal  flow  creates  an  interfacial  de¬ 
pression  of  the  pycnocline  and  generates  an  internal  baroclinic 
bore  at  the  sill  during  maximum  flows  in  a  tidal  cycle  ([12], 
[6],  [13],  [1],  [2],  [14].  In  the  Strait  of  Messina,  during  rema 
montante  when  the  heavier  LIW  crosses  the  sill,  the  pycno¬ 
cline  is  moved  upwards  at  the  sill  and  depressed  in  the  north 
basin.  The  depression  is  an  internal  bore  that  generates  a  south¬ 
wards  and  a  northwards  propagating  internal  wave.  The  north¬ 
wards  propagating  bore  steepens  and  disintegrates  into  solitary 
waves  whereas  the  southwards  propagating  bore  is  stopped  by 
the  sill.  When  the  semidiurnal  tide  reverses  to  rema  scendente^ 
the  southwards  propagating  internal  bore  of  depression  under¬ 
goes  a  hydraulic  jump  over  the  sill  and  into  the  southern  basin 
where  it  propagates  away  from  the  sill.  As  it  propagates,  non¬ 
linear  effects  steepen  its  leading  edge  until  it  disintegrates  into 
a  series,  or  “train,”  of  interfacial  nonlinear  short  internal  waves 
because  of  frequency  and  amplitude  dispersions.  The  process 
is  very  well  described  in  [2].  The  Aaracteristics  for  internal 
waves  in  the  Strait  of  Messina  are  a  propagation  speed  between 
0.80  and  1.00  m  s“^,  with  oscillations  up  to  2  amplitude 
(20-50  m).  When  the  train  of  internal  waves  is  well  formed, 
4-10  periods  ranging  from  8  to  30  rtiin  can  be  observed,  cov¬ 
ering  an  average  total  duration  of  2  h  [7]. 

B.  Presentation  of  the  Oceanographic  Data 

One  thermistor  chain  was  bottom-moored  about  27  km  south 
of  the  sill  of  strait  (Fig.  I ),  at  location  (1 5®  3  8' 5  7"  E,  38®02'24" 
N)  from  November  4, 2006  1 6:48:00  LJTC  to  November  7, 2006 
07:12:00  UTC,  with  a  1-min  temporal  resolution  which  allows 
clean  observations  of  internal  solitary  wave  events.  It  had  ten 
sensors  positioned  from  14-  to  128-m  depth. 

Fig.  2  shows  the  temperature  measured  by  the  thermistor 
string.  The  thermocline  depth  stands  around  95  m  for  the  first 
day  and  half,  and  raises  slowly  to  a  60-m  depth  after  November 
6,  2006.  The  deepening  of  warm  surface  waters  illustrates  the 
occurrence  of  internal  solitary  wave  train  events  every  half 
day.  After  November  6,  2006,  the  mean  pycnocline  depth  goes 
up  by  30  m,  most  likely  due  to  sortie  wind  event  relaxation 
(hypothesis,  no  meteorological  data). 

Table  1  summarizes  precisely  the  Occurrence  of  the  internal 
wave  trains  (IWTs)  observed  on  the  thermistor  string  data  and 
the  reversal  of  the  tide  from  rema  tnontante  (northwards)  to 
rema  scendente  (southwards).  It  shows  that  IWTs  are  visible 
6  h  after  the  tide  reversal  from  North  to  South.  If  considering 
the  hypothesis  made  by  [7]  that  the  internal  bores  are  gener¬ 
ated  at  slack  waters  when  the  flow  reverses  from  north  to  south 
(for  southwards  propagating  internal  solitary  waves),  then  the 
average  speed  of  the  IWTs  here  is  1 .55  m  s”^.  The  amplitudes 
of  the  internal  solitary  waves  range  between  40  and  50  m.  The 
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Fig.  2.  Thermistor  string  measured  temperature  from  Novxmbcr  4,  2006  1 6:48:00  UTC  to  November  7,  2006  07:12:00  UTC.  The  color  bar  displays  the  temper¬ 
ature  in  degrees  Celsius.  The  16  ®C  temperature  contour  is  plotted  to  underline  the  temperature  variability  due  to  occurrence  of  internal  waves. 


Table  I 

Internal  Wave  Trains  (1  WT)  Identified  on  the  Thermistor  String  Measurement.  As  the  Bore  Generated  at  “Rema  Montante'"  on  the  Northern 
Side  of  the  Strait  Undergoes  a  Hydraulic  Jump  Over  the  Sill  to  the  Southern  Basin  When  the  Tide  Reverses  Southwards  (Rema  Scendente), 
THE  Time  Difference  Between  the  SAR  Image  and  the  Tidal  Reversal  Is  Computed  in  Order  to  Get  the  Propagation  Speed  of  the  IWT  in 
THE  South  Basin.  The  First  Line  Indicates  the  Date  and  Time  of  the  IWT  Event,  and  the  Second  Line  Indicates  the  Time  of  Tidal 
Reversal  From  North  to  South  (Slack  Waters)  Previous  to  the  SAR  Image.  The  Last  Line  Shows  th$  Time  Difference 
Between  the  Previous  Tide  Reversal  and  the  Occurrence  of  the  IWT 


Inlemal  Wave  Trains 
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5:50 

6:01 

6:27 

6:01 

trains  have  two  to  four  internal  solitary  waves.  The  first  and 
second  ones  are  very  neat  in  the  data. 

C  Theory  of  the  EOF  Analysis 

The  EOF  decomposition,  also  called  principal  component 
analysis,  is  one  of  a  larger  class  of  inverse  techniques  and  is 
equivalent  to  a  data  reduction  method  [15].  It  is  one  of  the  most 
powerful  methods  for  reducing  the  variability  in  the  data  into 
a  few  patterns.  It  provides  a  compact  description  of  the  spatial 
and  temporal  variability  of  data  series  in  terms  of  orthogonal 
functions,  or  statistical  “modes,”  resulting  from  the  decomposi¬ 
tion  of  the  data  signal  in  eigenvalues  and  eigenvectors.  Each  set 
of  eigenvector,  called  EOF  mode,  represents  a  standing  oscilla¬ 
tion  pattern  [16].  The  eigenvalues  are  normalized  and  ordered 
according  to  their  size  to  indicate  what  EOF  mode  contains 
more  variance.  Therefore,  each  eigenvalue  is  associated  with  a 
pattern,  giving  a  measure  of  the  fraction  of  the  total  variance 
explained  by  each  EOF  mode.  The  projection  of  the  data  onto 
the  standing  EOF  modes,  called  principal  component  time 
series  (PC)  or  expansion  modes  of  the  EOFs,  gives  an  estimate 
of  how  the  pattern  oscillates  in  time.  Most  of  the  variance  is 
usually  in  the  first  few  orthogonal  functions.  However,  one 
must  be  careful  before  making  any  physical  assumption  on 


the  mechanism  of  a  phenomenon.  As  the  method  finds  a  set 
of  orthonormal  basis  vectors  that  maximizes  the  projection  of 
the  data  on  the  basis  vectors,  the  EOF  modes  are  uncorrelated 
over  space  and  the  expansion  coefficients  are  uncorrelated  over 
time.  Thus,  these  patterns  of  variability  are  primarily  “data” 
modes  and  not  necessarily  “physical’^  modes  [17]. 

In  this  paper,  the  EOF  analysis  was  applied  on  the  vertical 
to  the  temperature  anomaly  (the  work  was  done  on  the  density 
anomaly  in  [5]),  so  that  the  measured  temperature  could 

be  written  as 

M 

ip{z,t)  =  ‘ip{z)  +  '^ai{t)-(t>i{z)  (1) 

1  =  1 

where  t?(z)  is  the  time-averaged  vertical  temperature  profile, 
M  is  the  number  of  considered  EOfJ  ai{t)  is  the  amplitude  of 
the  ith  orthogonal  mode  at  time  i,  aind  <t>i{z)  is  the  ith  EOF. 
This  equation  says  that  the  time  variation  ofthe  dependent  scalar 
variable  rp{Zyt)  at  each  depth  results  from  the  linear  combina¬ 
tion  of  M  spatial  functions  (pi,  whos0  amplitudes  are  weighted 
by  M  time-dependent  coefficients  oci[t).  The  weights  a,(t)  tell 
how  the  spatial  modes  (pi  vary  with  time.  As  there  are  many  ter¬ 
minologies  for  EOF  analysis,  we  insist  on  the  fact  that  the  use  of 
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EOFI  EOF2 


Second  expansion  coefficient 


Fig.  3.  (Top)  Amplitude  of  the  first  two  EOFs  and  <f>2  for  the  thermistor  string  temperature.  The  first  mode  does  not  change  sign  over  the  water  column. 
(Bottom)  Corresponding  expansion  coefficients.  The  first  expansion  coefficient  clearly  represents  the  oscillation  of  the  thermocline  in  time  (compare  with  16  ®C 
temperature  contour  in  Fig.  2).  Both  coefficients  have  stronger  (absolute)  values  during  internal  wave  events. 


“EOFs”  refers  to  the  vertical  spatial  patterns  (l>i{z)  and  that  “ex¬ 
pansion/projection  coefficient”  refers  to  the  temporal  patterns 
or  weights  Qi(i)- 

Fig.  3(top)  illustrates  the  first  two  EOF  modes,  and 
Fig.  3(bottom)  shows  the  corresponding  principal  compo¬ 
nents.  With  the  inherent  efficiency  of  this  statistical  description 
a  very  few  empirical  modes  are  used  to  describe  the  funda¬ 
mental  variability  in  a  very  large  data  set;  the  two  first  modes 
for  our  data  set  represent  99.5%  of  the  variance.  All  EOFs 
are  very  clean  and  smooth,  representing  a  dynamical  pattern 
consisting  of  a  background  stratification  interacting  with  in¬ 
ternal  wave  occurrence.  The  amplitude  of  the  first  EOF  does 
not  change  sign  over  the  water  column  and  has  one  extremum. 
Therefore,  all  the  isotherms  oscillate  in  the  same  direction 
[18]  with  maximum  displacement  at  85-90  m.  The  second 
EOF  has  two  extrema,  as  a  consequence  of  the  orthogonality 
constraint,  and  changes  sign  close  to  the  first  EOF  maximum 
depth.  This  indicates  that  isotherms  are  moving  in  opposite 


directions  above  and  below  the  pycnocline  [18].  The  first  EOF 
can  be  interpreted  as  a  vertical  displacement  of  the  thermocline 
whereas  the  second  one  illustrates  the  change  in  slope  of  the 
temperature  gradient.  To  make  sure  the  physical  interpretation 
is  right  as  the  EOF  data  modes  are  not  necessarily  related  to 
physical  modes,  the  corresponding  first  two  projection  (ex¬ 
pansion)  coefficients,  illustrated  in  Fig.  3(bottom),  are  also 
plotted.  The  projection  coefficient  of  EOFI  is  negative  at  the 
beginning  and  turns  positive  after  November  6,  2006.  Looking 
at  the  thermocline  on  Fig.  2,  it  is  noticeable  that  the  projection 
coefficient  represents  exactly  the  variations  of  the  thermocline, 
being  negative  when  the  thermoclinC  is  deeper  than  its  mean 
depth  (EOFI  maximum)  and  positive  when  the  thermocline  is 
shallower.  The  two  periods  identified  in  Fig.  2  are  thus  clearly 
identified  in  the  first  expansion  coefficient  representing  the 
vertical  displacement  of  the  thermocline  depth. 

The  second  projection  coefficient  mostly  stays  positive  but 
the  amplitude  of  its  variations  is  bigger  (local  extrema)  during 
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Crescent  shape  expansion  coefficient  scatter  diagram 


Fig.  4.  Scalier  diagram  formed  by  plotting  ihc  second  expansion  coefficieni  versus  the  firsi  one.  The  scatter  ploi  has  a  significant  crescent  shape  illustraiing  thai 
despile  the  orthonormal  constraints  of  the  EOF  analysis,  the  firsi  I  wo  modes  are  dynamically  linked.  The  letters  correspond  to  the  expansion  coefficients  of  the 
profiles  rebuill  in  Fig.  5.  The  black  curve  represents  the  fourth-order  polynomial  fil  lo  the  cloud  of  points. 


an  internal  wave  event.  All  these  evidences  illustrate  that  the 
EOFs,  in  our  case,  can  be  linked  to  physical  modes. 

D.  Expansion  Coefficient  Scatter  Diagram 

As  most  of  the  variance  was  described  in  the  first  two  EOF 
modes,  the  scatter  diagram  of  the  first  two  expansion  coeffi¬ 
cients  was  plotted  to  see  if  they  could  be  related  to  an  internal 
wave  event  [5].  The  resulting  cloud  of  points  has  a  crescent 
shape  (Fig.  4)  meaning  that  the  two  coefficients  are  not  inde¬ 
pendent:  a2(<)  =  /(ai(<)). 

The  IW  signal  can  then  be  fully  described  as 

-  rp  +  +  a2{t)<j>-2{z)  (2) 

=  ^  +  ai<l>i  +  f{ai)<l>2.  (3) 

Two  regimes  can  be  distinguished  on  the  scatter  diagram.  In  the 
center  of  the  cloud  of  points  [|ai(i)|  -<  15  and  a2(^)  -<  0], 
where  the  cloud  of  points  is  thick  and  dense,  Qi{t)  takes  small 
values  meaning  a  weak  translation  of  the  thermocline  depth. 
a2(<)  is  negative,  and  the  corresponding  profiles  are  smooth  and 
typical  of  a  well-stratified  water  column  [5].  This  regime  rep¬ 
resents  a  “rest  regime”  with  no  internal  wave  events.  Looking 
now  at  the  crescent  “branches”  with  an  asymptotic  behavior 
[|ai(<)|  y  15  and  a2(0  y  0],  both  coefficients  increase  in  a 
linear  correlation.  This  regime  represents  typically  the  internal 
wave  dynamics:  when  the  thermocline  is  depressed  by  the  inter¬ 
facial  waves  [represented  by  ai(<)],  the  deviation  of  the  inter¬ 
face  generates  a  horizontal  circulation  with  opposite  horizontal 
currents  below  and  above  the  interface  of  the  thermocline  [gra¬ 
dient  variability  represented  by  a2{t)]‘  These  baroclinic  cur¬ 
rents  are  added  to  the  external  barotropic  current.  Once  the  1 WT 
has  passed,  this  baroclinic  structure  disappears  and  we  are  back 
in  the  first  described  regime. 


This  crescent  shape  scatter  diagram  was  confirmed  by  numer¬ 
ical  simulations  [5]  using  the  nonlinear  and  nonhydrostatic  2.5 
D  Lamb  model  [6]. 

The  crescent-shaped  distribution  tneans  that  the  two  EOF 
projection  coefficients  are  dynamically  linked.  When  this  re¬ 
sult  was  first  published  in  [5],  a  lot  ot  questions  arose  from  the 
reviewers  and  the  scientific  community.  EOFs  are  built  under 
some  orthogonal  constraints:  1)  the  EOF  modes  (f>iiz)  are  or¬ 
thogonal  in  space  meaning  that  there  is  no  correlation  between 
any  two  EOFs  (0*  •  <l>j  =  Sij)  and  2)  the  expansion  coefficient 
time  series  Qi{t)  is  orthogonal  in  time  meaning  that  there  is 
no  simultaneous  temporal  correlation  between  any  two  expan¬ 
sion  coefficients  (a^  •  aj  =  Sii).  Tliese  two  constraints  were 
checked  and  respected  for  the  first  two  EOF  modes  considered 
here:  <l>i  •  (f>2  =  0  and  ai  •  Q2  =  0.  The  scatter  diagram  of 
the  expansion  coefficients  having  a  specific  shape  is  not  in  con¬ 
tradiction  with  the  orthogonality  constraint;  it  only  highlights 
a  dynamical  interaction  between  these  two  modes  illustrating 
that  the  vertical  oscillation  of  the  thermocline  during  an  internal 
wave  event  is  related  to  the  change  in  slope  of  its  gradient. 

The  black  curve  in  Fig.  4  indicates  a  fourth-order  polynomial 
fit  P{qi)  of  the  scatter  diagram  and  thus  the  mathematical  rela¬ 
tion  between  the  first  two  expansion  (Coefficients.  This  property 
was  used  to  rebuild  the  profiles  from  the  scatter  plot  polynomial 
fit  and  the  two  vertical  EOF  modes. 

Two  other  examples  of  scatter  diagrams  are  displayed  in 
Fig.  5.  These  are  scatter  diagrams  which  were  rebuilt  from 
the  expansion  coefficients  of  themrtistor  string  temperature 
time  series  measured  close  to  the  Setubal  canyon  on  the 
western  coast  of  Portugal  during  the  2007  Random  Array  of 
Drifting  Acoustic  Receivers  (RADAR)  trial  jointly  organized 
by  CINTAL,  Algarve,  Portugal  and  Instituo  Hidrografico 
(IH),  Lisboa,  Portugal,  with  three  other  teams  from  NATO 
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Fig.  5.  Two  other  examples  of  scatter  plots.  The  temperature  thermistor  strings  measurements  were  made  during  the  RADAR*07  trial  close  to  the  Setubal  canyon 
on  the  western  coast  of  Portugal  in  July  2007.  The  RADAR’07  trial  was  jointly  organized  by  CINTAL/SiPLAB  and  Instituo  H(drografico  (IH)  with  three  other 
teams  from  NURC,  ULS  Research,  and  Navy  Research  Laboratory  (NRL).  Their  positions  arc  (8®  58^6"  W;  38®  1 9' 44"  N)  for  thermistor  string  I .  and  (8®  5507" 
W;  38®2lM5"  N)  for  thermistor  string  2.  The  emettor  depths  of  hermistor  string  I  ranged  from  15  to  99  m  with  a  spacing  of  1  m.  Thermistor  string  2  had  1 1 
emettors  ranging  from  2  to  22  m  by  step  of  2  m.  They  were  both  positioned  northward  of  the  Setubal  canyon. 


Fig.  6.  Randomly  rebuilt  temperature  profiles  from  the  two  first  modes  of  the  EOF  analysis.  In  black,  original  data  temperature  profiles,  in  gray,  two-mode  rebuilt 
profiles.  The  profiles  are  located  on  the  expansion  coefficient  scatter  diagram  in  Fig.  4. 


Undersea  Research  Centre  (NURC),  La  Spezia,  Italy,  Navy 
Research  Laboratory  (NRL),  Stennis  Space  Center,  MS,  and 
Heat,  Light  and  Sound  Research  (HLS),  La  Jolla,  CA.  Their 
positions  are  (8°58'6"  W;  38°  1 9' 44"  N)  for  thermistor  string 
1,  and  (8°55'37"  W;  38°2r45"  N)  for  thermistor  string  2. 
The  internal  wave  dynamics  here  is  not  as  strong  as  in  the 
Messina  Strait  (weaker  amplitudes)  and  the  dynamics  close 
to  the  canyon  are  much  more  complex.  The  shape  of  the 
cloud  of  points  in  both  scatter  diagrams  is  not  as  neat  as  for 
Messina  Strait.  Still  the  lower  limit  of  the  cloud  of  points  is 
well  designed  and  can,  once  again,  be  fitted  with  a  polynomial 
fit  illustrating  a  dynamical  link  between  both  coefficients.  This 
study  will  go  on  with  the  Messina  Strait  scatter  plot,  but  these 
latter  examples  illustrate  that  our  methodology  can  be  derived 
for  several  areas  with  strong  dynamics. 


Ill.  Representation  of  Temperature  With  EOFs 

In  [5],  the  EOF  analysis  was  done  on  the  density.  As  our  aim 
here  is  to  rebuild  the  temperature  profiles  for  sound-speed  pro¬ 
file  computation,  the  EOFs  are  applied  to  the  temperature  mea¬ 
sured  by  the  thermistor  string.  The  EOF  and  the  scatter  diagram 
are  obviously  exactly  the  same;  only  the  axis  scales  change. 


A.  Reconstruction  Methodology 

As  noticed  in  Section  II-C,  the  variance  is  fully  contained  in 
the  first  two  EOFs  (over  99.5%).  Moreover,  the  EOF  coefficients 
are  dynamically  linked  and  the  whole  water  column  variability 
can  be  described  by  a  polynomial  function. 
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Crescent  shape  expansion  coefficient  scatter  diagram 
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Fig.  7.  (Top)  Scatter  diagram  formed  by  plotting  the  second  expansion  coefficient  versus  the  first  one.  The  black  dots  represent  the  sampling  in  20  profiles  of  the 
polynomial  fit.  The  corresponding  rebuilt  profiles  are  illustrated  below,  from  the  first  step  of  the  two-mode  reconstruction  to  the  computation  of  sound  speed. 


•  By  choosing  a  point  on  the  polynomial  fit  P{ai)  of  the 
scatter  diagram,  the  water  column  temperature  can  be  re¬ 
constructed  using 

Tiiz)  —  ^mean  +  an  *  4>i{z)  +  P(aii)  *  4>2(z), 

for  i  =  time  (4) 

where  4>i  and  <(>2  are  the  EOFs  (spatial  modes)  and  qi^  and 
P{aii)  are  the  coefficients  at  the  chosen  time  i. 


Moreover,  as  the  thermistor  string  data  were  available  between 
20-  and  1 30-m  depth,  the  temperature  profiles  were  then  extrap¬ 
olated  on  the  whole  water  column  (st^s  in  Fig.  7). 

•  For  the  layer  up  to  the  surface,  the  profile  was  extended 
by  a  constant  value  corresponding  to  the  shallower  value 
available.  As  we  are  in  the  mixed  layer,  this  approximation 
is  correct. 

•  From  the  deepest  profile  value  to  the  bottom,  the  profile 
was  extrapolated  to  a  mean  temperature  profile  computed 
from  the  deepest  CTD  measurements.  The  profiles  were 
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then  slightly  filtered  to  smooth  the  transitions  of  the  dif¬ 
ferent  parts  of  each  profile  (Fig.  7,  middle  graph). 

Four  profiles  were  rebuilt  randomly  and  are  indicated  in  Fig.  4. 
Profiles  A  and  B  are  in  the  “rest  regime”  where  the  density  of 
points  is  most  dense.  Profiles  C  and  D  are  located  farther  along 
the  asymptotes,  corresponding  more  to  an  internal  wave  regime. 
For  profiles  located  in  the  rest  regime  (first  two  graphs  in  Fig.  6), 
the  reconstruction  works  well.  A  slight  temperature  inversion 
is  observed,  most  likely  due  to  the  very  negative  values  of  the 
second  mode  expansion  coefficient.  Unfortunately  for  the  last 
reconstructed  profile  [Fig.  6(right)],  the  result  is  poor.  To  com¬ 
pensate  for  the  high  therm ocline  depth  (60  m),  the  expansion  co¬ 
efficients  are  strong,  thus  introducing  too  much  variability  and 
a  great  artefact  at  the  theimocline  depth  mean  position  ( 1 1 0  m). 
This  was  corrected  by  a  temperature  inversion  check  in  the  re¬ 
construction  routine. 

•  The  temperature  inversion  check  is  simple.  If  the  tempera¬ 
ture  value  at  level  Zi^i  (deeper)  is  higher  than  the  temper¬ 
ature  at  level  Zi  then  the  temperature  at  Zi^i  is  corrected 
and  taken  equal  to  the  value  at  That  simple  correction 
works  well  as  the  inversions  often  appeared  in  the  mixed 
layer  or  at  the  bottom  of  the  theimocline. 

•  Finally,  the  sound  speed  c^(^)  was  computed  using 
Medwin  formula  from  [19],  function  of  temperature, 
salinity  (salinity  was  chosen  from  data  as  constant  and 
equal  to  38.51 ),  and  depth 

Ci{z)  =  1449.2  +  4.6  x  Ti{z)  -  0.055  x  T^{z)  +  ■■■ 
0.00029  X  T^iz)  +  (1.34  -  Ti(z))  x  {S  -  35) 

+  0.016X2.  (5) 

These  steps  are  illustrated  in  Fig.  7.  The  sound-speed  profiles 
are  now  ready  to  be  used  as  acoustical  model  inputs. 

B.  Case  I:  Curve  Sampling  and  Systematic 
Reconstructed  Profiles 

By  sampling  the  polynomial  fit,  an  exhaustive  set  of  temper¬ 
ature  profiles  can  be  reconstructed.  The  curve  was  sampled  for 
20  values  of  the  first  principal  component.  It  was  chosen  not 
to  start  at  the  minimum  of  the  first  expansion  coefficient  as  the 
density  of  points  was  very  low  and  most  likely  insufficient  to 
have  the  proper  shape  of  the  polynomial  left  branch  (Fig.  7). 

Once  the  sampling  coefficients  were  chosen,  the  20  sound- 
speed  profiles  were  reconstructed  using  the  methodology  de¬ 
scribed  in  the  previous  paragraph.  In  50%  of  the  cases,  the  fit 
to  the  data  is  very  good.  The  main  differences  appear  for  ex¬ 
treme  values  of  ai  (most  likely  because  of  the  weaker  density 
of  points  that  may  have  an  impact  on  the  polynomial  fit)  but  also 
when  the  gradient  of  the  theimocline  is  very  strong  as  in  profiles 
10  and  13.  Overall,  the  fit  is  good  for  a  two-EOF  reconstructed 
profile. 

These  profiles  illustrate  the  whole  water  column  variability 
that  can  be  observed.  Applying  an  acoustic  propagation  model 
to  these  profiles  will  give  an  exhaustive  representation  of  the 
environmental  acoustic  properties  of  the  Strait  (Section  IV) 
without  any  central  processing  unit  (CPU)  consuming  ocean 
modeling. 


CrMC«nl  slMp#  •xp«nsion  coefficient  scatter  diagram 
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Fig.  8.  Reconstruction  of  an  IWT.  The  color  bar  indicates  the  temperature  in 
degrees  Celsius.  The  profiles  were  chosen  on  the  polynomial  fit  in  order  lo  re¬ 
build  an  IWT  whose  characteristics  corresponded  to  the  ones  identified  in  [5J. 


C.  Case  2:  IWT  Reconstruction 


In  this  second  case,  we  have  looked  at  how  the  polynomial 
fit  was  browsed  during  an  IWT  event  A  wave  train  was  recon¬ 
structed  in  agreement  with  the  IWT  characteristics  in  the  Strait 
of  Messina. 

Fig.  8  illustrates  all  the  steps  from  the  sampling  of  the  profiles 
to  the  reconstruction  of  the  IWT.  The  rebuilt  wave  train  has  four 
internal  waves  of  maximum  amplitude  of  30  m  in  a  distance  of 
4  km.  It  was  rebuilt  with  1 4  profiles. 

The  profile  at  the  sill  is  the  one  that  has  the  shallower  thermo- 
cline,  so  it  was  chosen  well  on  the  “right-hand  side”  (positive 
first  expansion  coefficient)  of  the  polynomial.  The  profile  south¬ 
wards  must  have  a  deeper  thermoclint  and  was  chosen  closer  to 
the  bottom  of  the  polynomial  to  have  k  theimocline  about  1 00  m 
deep.  Concerning  now  the  train  itself,  points  2,  4,  and  6  repre¬ 
sent  the  main  waves  of  the  train,  and  we  must  then  have  negative 
values  of  the  first  expansion  coefficient  to  get  a  deeper  thermo- 
cline:  the  values  of  the  first  expansion  coefficient  and  thus  the 
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Fig.  9.  Two-dimensional  transmission  loss  on  the  whole  water  column  for  /  =  1 00  Hz  (top)  and  /  =  1 500  Hz  (bottom)  for  Profile  8  “rest-position*’  (left  column) 
and  Profile  1  “internal  wave  occurrence”  (right  column).  Color  bars  indicate  the  loss  in  decibels  for  each  configuration.  The  range  goes  from  the  Sicilian  coast 
to  the  Calabrian  coast.  On  the  first  and  third  lines,  the  source  is  20  m  deep,  and  on  the  second  and  fourth  lines,  the  source  is  1 30  tn  deep.  The  range-independent 
sound-speed  profile  is  displayed  in  the  graphs  in  black. 


amplitude  of  each  wave  are,  respectively,  20,  1 7,  and  10  m.  Be¬ 
tween  each  wave,  the  thermocline  has  to  come  back  to  its  av¬ 
erage  depth,  so  that  the  profiles  are  chosen  close  to  the  bottom 
of  the  polynomial. 

Just  by  using  the  polynomial  function,  a  wave  train  can  be 
easily  reconstructed  and  thus  a  physical  feature  study  can  be 
conducted. 

IV.  Acoustic  Propagation 

In  this  section,  we  illustrate  the  methodology  by  running 
acoustic  simulation  across  and  along  the  straits  with  recon¬ 


structed  environments.  The  purpose  is  not  to  derive  a  full 
acoustic  study  of  propagation  in  the  Strait  of  Messina  but  rather 
to  illustrate  by  examples  the  potentiality  of  our  reconstruction 
approach.  The  acoustical  propagation  runs  are  obtained  with 
the  range-dependent  acoustic  model  (RAM),  developed  in  [20]. 
This  model,  based  on  split-step  Pad^  resolution  of  the  parabolic 
equation  (PE),  allows  the  use  of  a  wide  range  of  frequencies.  It 
is  used  hereafter  from  50  to  3500  Hz  ^y  a  step  of  200  Hz. 

Two  configurations  were  chosen i  one  was  implemented 
across  the  strait  (transverse  configufation  with  sound-speed 
profiles  described  in  Case  I  in  Section  lll-B)  and  the  other 


754 


IEEE  JOURNAL  OF  OCEANIC  ENGINEERlNd.  VOL.  36.  NO.  4,  OCTOBER  2011 


DOWNSLOPE  20M  <  Frequency  x  1 0OHz  DOWNSLOPE  20M  *  FrtquwKy  x  2500Hz 


DOWNSLOPE  130M  -  Fr*qtMncy  x  lOOHz 


Sourc*  /  R«c«lv«r 
Dvplhs  X  I30ni 


DOWNSLOPE  130U  <  Fr^quvncy  -  2500Hz 


30 

Rang*  (km) 


—  WlthlW 

—  HotW 


23  25  27  30 

Rang*  (kfn) 


Fig.  10.  Transmission  loss  between  a  source  and  a  receiver  in  the  axial  downslope  configuration  (along  the  median  of  the  Strait)  for  frequencies  100  Uz  (left 
column)  and  2500  Hz  (right  column).  Source  is  positioned  9  km  away  from  the  sill  (in  front  of  Messina  town),  and  loss  is  dispteyed  from  20  to  40  km  (so  that 
the  sill  is  on  the  left  of  the  graph  and  Ionian  Sea  on  the  right).  The  IWT  goes  from  23  to  27  km.  In  the  first  line,  the  source  and  the  receptor  are  at  the  20-m  depth 
(study  of  the  surface  duct).  In  the  second  line,  the  source  and  the  receptor  are  at  the  130-m  depth  (below  the  thermocline  and  the  $ound-spced  minimum). 


one  followed  the  median  of  the  strait  (axial  configuration  with 
sound-speed  field  described  in  Case  2  in  Section  III-C). 

•  In  the  configuration  transverse  to  the  strait,  the  source  is  po¬ 
sitioned  900  m  away  from  the  Sicilian  coast  and  emits  to¬ 
wards  the  Calabrian  coast  ( 1 6  600-m  distance).  The  scatter 
plot  polynomial  is  sampled  in  20  points  and  runs  were 
made  systematically  for  every  reconstructed  profile  and 
every  frequency:  the  sound-speed  profile  is  range  indepen¬ 
dent  but  not  the  sediment  nature  (see  below). 

•  In  the  axial  configuration  along  the  strait  median,  the 
source  is  positioned  either  in  front  of  Messina  (9  km  south 
of  the  sill)  and  emits  downslope  until  40  km  south,  or  40 
km  south  of  the  sill  and  emits  upslope  towards  the  sill. 
The  sound-speed  field  contains  a  reconstructed  IWT  and 
is  thus  range  dependent. 

Bottom  nature  is  derived  from  [2 1  ]  and  from  data  extracted  from 
SHOM  sediment  database.  For  the  transverse  configuration,  the 
bottom  nature  is  a  succession  of  rocks  (from  the  source  to  6  km), 
coarse  sand  and  silt  (until  5  km),  silt  (until  1 1  km),  and  finally 
coarse  sand  and  silt  (until  the  receiver).  For  the  axial  configura¬ 
tion,  the  bottom  consists  in  a  section  of  coarse  sand  and  silt  over 
the  first  9  km  and  silt  elsewhere. 

The  source  depths  are,  respectively,  20  m  (shallow,  in  the 
surface  duct),  90  m  (intermediate,  in  the  thermocline),  and  130 
m  (deep,  in  the  channel). 

A.  Transverse  Configuration 

The  transverse  configuration  is  based  on  the  Case  1  profile 
construction.  The  acoustical  model  is  positioned  transverse  to 
the  strait,  27  km  south  of  the  sill,  and  the  sound-speed  profile  is 
range  independent.  The  used  profiles  are  the  20  profiles  rebuilt 


as  explained  in  Section  III-B.  For  each  of  these  profiles,  propa¬ 
gation  loss  simulations  were  made  for  frequencies  ranging  from 
50  to  800  Hz  by  a  step  of  50  Hz  and  from  900  to  3500  Hz  by  a 
step  of  200  Hz. 

An  example  of  transmission  loss  field  is  displayed  in  Fig.  9. 
Profile  8  is  located  at  the  bottom  of  the  scatter  diagram  (see  po¬ 
sition  and  sound-speed  profile  in  Fig.  7),  thus  corresponding  to 
a  “no  internal  wave/rest  position”  profile  (left  column).  Profile 
I  represents  the  occurrence  of  an  internal  wave  event  deepening 
the  thermocline  (right  column).  The  frequencies  displayed  are 
1 00  (four  top  pictures)  and  1 500  Hz  (four  bottom  pictures).  The 
sound-speed  profiles  are  displayed  on  each  graph.  First  and  third 
lines  represent  the  transmission  loss  for  a  20-m-deep  source, 
and  second  and  fourth  lines  represent  the  transmission  loss  for 
a  130-m-deep  source. 

Several  well-known  conclusions  easily  come  up. 

•  Sediment:  for  the  lower  frequencies  (and  thus  longer  wave¬ 
lengths),  the  nature  of  the  sediment  can  induce  bottom  loss. 
Whereas  this  loss  is  little  for  rocky  bottoms,  the  acoustic 
energy  penetrates  deeper  for  softer  sediments.  For  higher 
frequencies  (bottom  pictures),  the  signal  does  not  penetrate 
at  all  in  the  sediment.  It  also  affects  the  interference  pat¬ 
terns  at  low  frequencies.  The  rays  that  have  first  touched 
the  rocky  bottom  have  a  much  finer  interference  pattern 
than  the  rays  reflecting  on  the  silt  bottom  or  being  guided 
in  the  water  column.  For  the  high  frequencies,  the  change 
in  sediment  nature  has  hardly  apy  impact. 

•  Surface  duct:  for  both  sound-speed  profiles,  the  gradient 
close  to  the  surface  is  slightly  positive  because  of  the 
isothermal  mixed  layer.  This  bit  is  followed  by  a  strong 
negative  gradient.  When  the  source  is  positioned  in  the 
positive  gradient  (20-m  depth),  the  curvature  of  the  rays 
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Fig.  II.  Two-dimensional  transmission  loss  illustrating  the  case  from  Fig.  10  when  the  source  and  the  receiver  are  at  the  l30-m  depth,  with  (top)  or  without 
(bottom)  IWT.  The  left  column  displays  the  results  for  the  first  300  m  of  the  water  column  at  frequency  100  Hz,  and  the  right  column  at  frequency  2500  Hz.  Color 
bars  indicate  the  transmission  loss  in  decibels. 


is  upward.  The  rays  with  a  low  angle  propagate  in  this 
channel  by  multiple  reflections  on  the  surface.  The  trans¬ 
mission  loss  is  usually  quite  low  if  the  frequency  is  neither 
too  low  (low-frequency  cutolT)  nor  too  high  (surface  re¬ 
flection  loss).  For  the  1 00-Hz  run,  the  energy  leaks  out  of 
the  duct  after  approximately  4  km.  For  the  1500-Hz  run, 
the  energy  is  guided  through  the  whole  strait. 

•  Shadow  zone:  for  the  20-m  depth  source,  the  rays  with 
a  higher  angle  than  the  ones  propagating  in  the  surface 
duct  are  quickly  going  out  of  the  duct  and  their  curvature 
is  strongly  deviated  to  the  bottom  because  of  the  strong 
negative  sound-speed  gradient.  At  1500  Hz,  this  leaves  an 
important  shadow  zone  which  can  only  be  and  is  insonified 
by  bottom  reflected  rays  (which  is  less  efficient  than  direct 
rays).  At  1 00  Hz,  the  energy  leaking  out  of  the  surface  duct 
is  refracted  on  the  deeper  sound-speed  positive  gradient, 
which  allows  a  more  important  insonifleation  of  the  water 
column  than  observed  at  1 500  Hz.  A  reinsonification  of  the 
surface  duct  is  thus  observed  at  12  km. 

•  Deeper  source:  when  the  source  is  placed  at  130  m,  the 
surface  duct  is  not  insonified  anymore  except  by  some 
geometrically  reflected  rays.  For  profile  8,  the  source 
immersion  is  close  to  the  sound-speed  minimum.  This 
corresponds  to  the  best  insonifleation  possible  [SOund 
Frequency  And  Ranging  channel  (SOFAR)  propagation 
type — the  SOFAR  channel  acts  as  a  waveguide  for  sound 


and  sound  waves  within  the  channel  may  travel  long 
distances  before  dissipating]  thanks  to  the  absence  of 
reflection  loss  and  the  presence  of  multiple  tracks  (the 
only  loss  is  water  absorption).  It  is  well  illustrated  at 
1500  Hz.  The  strait  width  is  too  narrow  to  observe  any 
convergence  zone.  For  profile  1,  the  source  position  is 
under  the  isothermal  zone  but  still  in  the  thermocline 
(above  the  sound-speed  minimum).  In  that  case,  a  huge 
shadow  zone  is  observed  as  all  the  acoustical  energy  is 
deviated  downward. 

This  arbitrary  example  illustrates  the  different  phenomena  (sur¬ 
face  duct,  sofar  propagation,  sediment  impact)  that  can  take 
place.  The  results,  highly  depending  on  what  profile  is  used, 
highlight  the  crucial  importance  of  taking  into  account  internal 
wave  dynamics. 

B.  Axial  Configuration 

The  reconstructed  IWT  feature  i$  studied  with  either  the 
source  close  to  the  sill  of  the  strait  tind  emitting  towards  the 
open  sea  (downslope  propagation)  of  the  source  farther  down 
the  strait  and  emitting  towards  the  sill  (upslope  propagation). 

1)  Downslope  Propagation:  In  this  configuration,  we  look  at 
the  impact  of  a  complete  wave  train  along  the  axe  of  the  strait. 

The  source  is  positioned  in  front  of  Messina  (about  9  km 
south  of  the  sill);  the  receptor  is  positioned  downslope  after  the 
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Fig.  1 2.  Transmission  loss  between  a  source  and  a  receiver  in  the  axial  upslope  configuration  (along  the  median  of  the  Strait)  for  frequencies  1 00  Hz  (left  column) 
and  2500  (right  column).  Source  is  positioned  40  km  away  from  the  sill  (in  the  Ionian  Sea)  and  emits  towards  the  sill  thus  facing  the  IWT.  Loss  is  displayed  from 
—40  to  —20  km  (so  that  the  sill  is  on  the  right  of  the  graph  and  Ionian  Sea  on  the  left;  contrary  to  Fig.  1 0).  The  IWT  goes  from  ^27  to  —23  km.  In  the  first  line, 
the  source  and  the  receptor  arc  at  the  20-m  depth  (study  of  the  surface  duct).  In  the  second  line,  the  source  and  the  receptor  are  at  the  I30-m  depth  (below  the 
thcrmoclinc  and  the  sound-speed  minimum). 


wave  train  (40  km  away  from  the  strait).  The  internal  wave  event 
takes  place  from  23  to  27  km. 

Transmission  loss  is  displayed  versus  range  in  Fig.  10.  At 
low  frequencies  (left  column  in  Fig.  1 0),  there  is  hardly  any  dif¬ 
ference  between  the  “internal  wave”  case  and  the  “no  internal 
wave”  case.  So  the  leakage  observed  in  the  2-D  graphs  is  inde¬ 
pendent  from  the  internal  wave  presence.  Acoustic  propagation 
at  low  frequencies  is  dominated  by  the  sediment  nature  rather 
than  the  water  column  characteristics.  At  high  frequencies  (right 
column  in  Fig.  10),  the  IWT  impact  (extending  from  23  to  27 
km)  is  visible  right  after  the  event  (from  27  to  33  km)  when  the 
source  is  at  the  depth  of  the  event  or  deeper  (90,  not  shown,  and 
130  m). 

In  general,  the  presence  of  an  internal  wave  event  increases 
the  transmission  loss,  up  to  10  dB  in  the  case  of  surface  duct 
configuration  or  when  the  source  is  positioned  at  the  same  depth 
or  deeper  than  the  internal  wave  event  for  a  receptor  under  the 
mixed  layer  (130  m). 

To  understand  better  the  increased  transmission  loss  in  the 
presence  of  internal  waves  for  higher  frequencies,  Fig.  1 1  dis¬ 
plays  the  2-D  transmission  loss  at  1 00  and  2500  Hz  for  a  source 
130  m  deep  for  “internal  wave”  case  (bottom  line)  and  “no  in¬ 
ternal  wave”  case  (top  line).  Looking  at  the  2500-Hz  simulations 
(right  column  of  Fig.  1 1 ),  the  IWT  is  followed  by  an  energy  loss 
at  30  km,  which  explains  the  increase  of  transmission  loss  in 
Fig.  10.  As  the  wavelength  is  very  short,  the  acoustic  signal  sees 
the  train  as  one  thermocline  depression.  The  rays  are  refracted 
towards  the  bottom  thus  leading  to  a  shadow  zone  going  from 
the  first  solitons  at  27-33  km,  when  the  first  refracted  rays  get 
back  to  the  thermocline.  Without  this  event,  the  water  column 
is  fully  insonified  as  Fig.  1 1  (top  right)  shows. 

2)  Upslope  Propagation:  The  configuration  is  the  same  as 
the  previous  axial  one,  but  the  source  is  positioned  at  the  end 


of  the  strait  and  emits  towards  the  sill  so  that  the  propagation 
occurs  upslope. 

Fig.  12  (corresponding  to  the  downslope  configuration 
Fig.  1 0)  shows  the  same  conclusions  iiS  for  the  downslope  case, 
but  with  much  less  difference  between  the  simulations  with  and 
without  internal  waves.  Fig.  13  illustrates  the  upslope  acoustic 
propagation  for  a  130-m-deep  source  and  two  frequencies: 
100  Hz  (left  column)  and  2500  Hz  (right  column).  As  already 
noticed,  there  is  little  impact  on  low  frequencies.  On  the  con¬ 
trary,  at  2500  Hz,  the  deviation  of  th0  direct  ray  by  the  leading 
wave  is  clearly  visible.  It  is  followed  by  a  small  shadow  zone 
before  the  first  resurgence.  After  the  wave  train  the  direct  ray 
has  disappeared,  but  as  the  propagation  occurs  upslope,  the 
refracted  and  bathymetry-reflected  rays  do  insonify  the  water 
column  until  the  strait  and  the  differences  are  not  as  dramatic 
as  they  were  in  downslope  configuration. 

C.  Discussion 

The  EOF-reconstruction  method  enables  to  build  data-driven 
synthetic  environment  which  can  be  used  to  model  the  acoustic 
propagation  as  illustrated  previously.  This  enables  to  con¬ 
duct  parametric  studies  or  statistical  $tudies  by  multiplying  the 
number  of  simulated  environments  and  adapting  the  operational 
scenarios  to  the  addressed  problems  (e.g.,  ASW,  tomography 
experiment  design,  acoustic  communications  optimization, 
marine  mammals  risk  mitigation). 

As  an  illustration,  Figs.  14  and  15  give  an  example  of  poten¬ 
tial  tactical  aid  by  covering  the  whole  frequency  range  summing 
the  differences  between  the  simulations  with  and  without  in¬ 
ternal  waves  at  each  frequency  in  the  axial  configuration.  Fig.  1 4 
shows  that  in  the  upslope  case,  the  acoustic  propagation  is  de¬ 
viated  by  the  envelop  of  the  IWT  sO  that  a  quiet  area  is  cre¬ 
ated  below  the  thermocline  (green  area).  On  the  opposite  side. 
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Fig.  1 3,  Two-dimensional  transmission  loss  illustrating  the  case  from  Fig.  1 2  when  the  source  and  the  receiver  are  at  the  IBO^m  depth,  with  (top)  or  without 
(bottom)  I WT.  The  left  column  displays  the  results  for  the  first  1 000  m  of  the  water  column  at  frequency  100  FIz  (to  give  the  exandple  of  the  bottom  reflections  in 
the  upslope  configuration),  and  the  right  column  at  frequency  2500  Hz.  Color  bars  indicate  the  transmission  loss  in  decibels. 
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Fig.  14  Upslope  qualitative  insonification  map.  The  map  was  obtained  by 
adding  all  the  transmission  loss  differences  between  the  simulations  with  in¬ 
ternal  waves  and  without,  for  all  frequencies  from  50  to  2500  Hz.  The  green 
(resp.,  red)  color  means  the  transmission  loss  with  internal  wave  presence  is 
higher  (resp.,  lower)  than  without  IW. 


Other  areas  (in  red)  are  more  likely  to  be  insonified  in  case  of  an 
IWT.  Fig.  15  shows  that,  for  the  downslope  case,  low  frequen¬ 
cies  do  not  show  any  impact  of  the  water  column  physics,  but 
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Fig.  15.  Transmission  loss  versus  frequency  fbr  the  axial  downslope  configu¬ 
ration  for  a  1 30-m-deep  source.  The  transmission  loss  is  averaged  on  the  area 
defined  from  28  to  31  km  in  range  and  1()()-15D  m  in  depth)  in  the  presence  of 
internal  waves  (black  line)  and  without  internal  waves  (gray  line). 


from  400  Hz  and  on,  the  difference  becomes  important  because 
of  the  shadow  zone  induced  by  the  deflection  of  the  rays  by  the 
IWT  envelop. 

These  “statistical”  outputs  summar^e  the  impact  of  the  water 
column  on  any  frequencies.  They  are  CPU  cheap,  quickly  avail- 
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able,  and  can  be  produced  for  any  scenario,  which  gives  an  in¬ 
teresting  strategic  approach  to  this  methodology. 

V.  Conclusion 

This  paper  has  presented  a  novel  approach  to  synthesize  re¬ 
alistic  environments  for  ocean-acoustic  parametric  studies.  An 
EOF  decomposition  is  applied  to  a  temporal  series  of  temj^er- 
ature  profiles  on  which  it  can  be  observed  that  the  first  two 
time-dependent  expansion  coefficients  are  dynamically  linked 
by  the  internal  wave  effects.  The  expansion  coefficients  plus 
the  polynomial  can  be  used  to  reconstruct  temperature  profiles 
independently  from  the  set  of  data  and  can  be  used  to  infer 
the  acoustic  environment  in  the  area.  The  methodology  pre¬ 
sented  in  this  paper  is  simple  to  run  and  requires  a  very  afford¬ 
able  set  of  data.  It  could  be  used  as  an  efficient  alternative  to 
ocean  and  acoustic  model  coupling  for  process  studies  or  for  re¬ 
gional  studies  especially  in  poorly  known  areas  or  highly  vari¬ 
able  areas,  where  it  is  difficult  to  obtain  good  sound-speed  pro¬ 
file  prediction  from  ocean  models. 

This  method  could  be  of  great  interest  for  the  development  of 
statistical  ocean-acoustic  decision  aids  for  underwater  acoustic 
of)erations  (e.g.,  ASW,  tomography,  communications).  Initially 
a  set  of  data  of  an  area  is  needed.  The  statistical  analysis  gives 
a  polynomial  characterization  of  the  dynamics  in  the  area.  This 
could  be  derived  for  many  areas  with  strong  dynamics.  Once 
this  fX)lynomial  characterization  is  established,  the  users  do  not 
need  the  data  anymore  (gain  of  time  and  memory  storage)  and 
can  rebuild  any  possible  profile/scenario  for  physical  analysis 
or  statistical  tools,  possibly  through  a  “blackbox”  doing  the  re¬ 
building  work.  An  acoustic  propagation  model  can  then  be  cou¬ 
pled  to  this  black  box  to  simulate  any  desired  scenario.  Thanks 
to  this  method,  the  simulation  of  scenarios  and  the  choice  of  de¬ 
cision  aids  are  in  the  power  of  the  user  who  then  decides  exactly 
what  he  needs  according  to  one  situation. 

The  results  in  the  Strait  of  Messina  prove  that  statistical  char¬ 
acterizations  of  ocean  processes  could  be  an  interesting  tool  to 
replace  or  more  likely  complement  3-D  ocean-acoustic  predic¬ 
tions.  Future  work  will  investigate  the  possibility  to  extend  the 
EOF  reconstruction  method  to  other  oceanographic  phenomena 
such  as  fronts  and/or  eddies. 
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